Dynamic polarizability of rotating particles in electrorheological fluids 
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A rotating particle in electrorheological (ER) fluid leads to a displacement of its polarization 
charges on the surface which relax towards the external applied field Eo, resulting in a steady- 
state polarization at an angle with respect to Eo. This dynamic effect has shown to affect the 
ER fluids properties dramatically. In this paper, we develop a dynamic effective medium theory 
(EMT) for a system containing rotating particles of finite volume fraction. This is a generalization 
of established EMT to account for the interactions between many rotating particles. While the 
theory is valid for three dimensions, the results in a special two dimensional configuration show that 
the system exhibits an off-diagonal polarization response, in addition to a diagonal polarization 
response, which resembles the classic Hall effect. The diagonal response monotonically decreases 
with an increasing rotational speed, whereas the off-diagonal response exhibits a maximum at a 
reduced rotational angular velocity ooo comparing to the case of isolated rotating particles. This 
implies a way of measurement on the interacting relaxation time. The dependencies of the diagonal 
and off-diagonal responses on various factors, such as cjo, the volume fraction, and the dielectric 
contrast, are discussed. 

Keywords: Electrorheological fluid; Dynamic polarizability; Effective medium theory; Off-diagonal response 



INTRODUCTION 



In many aspects, an electrorheological (ER) fluid con- 
taining suspensions at rest is quite different from an ER 
fluid containing suspensions subjected to rotational mo- 
tions, [1, 2, 3, 4, 5, 6, 7, 8] which has been referred to 
the dynamic effect in ER fluid, [9, 10, 11] and is of great 
relevance in various applications of ER fluid. Reports 
regarding this dynamic effect have been extensive in the 
literature. [12, 13, 14, 15, 16, 17, 18] Typically there are 
three cases [11] for rotation of the suspending particles 
about their centers: (1) rotation due to external torque, 
which, for instance, might be induced by the shear flow 
in the ER fluid. [3, 9, 10, 11, 12, 13] (2) particle rota- 
tion due to a rotating external field, which is a common 
phenomena utilized in electrorotation assay. [14, 15] (3) 
spontaneous rotation of particular particles within weak 
conducting fluid in a dc field above a threshold . This 
is first discovered almost a century ago, [19] now recog- 
nized as Quincke effect or negative ER effect. [16, 17] In 
all these cases, the rotation axis is favored in a special 
direction. 

The rotational motion leads to a displacement of sur- 
face charge (not for Quincke rotation) on the particles 
that mediates the interactions in the ER system. For 
an isolated rotating particle, it has been shown that the 
displaced surface charge on the particle due to rotation 
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will relax towards the external applied field Eo- [9] The 
result of competition between the displacement and the 
relaxation is a steady-state polarization which deviates 
from the applied field's direction, and depends on the 
product of the rotational angular velocity ujq of the par- 
ticle and the relaxation time tq. At the same time, the 
magnitude of the steady-state polarization is reduced by 
a factor as compared to the resting case. [5, 9] which 
is determined by the same competition. As the inter- 
particle force depends crucially on the magnitude and 
direction of polarizations, the complications due to the 
rotational motion were shown to be dramatic. [10, 18] 
A first attempt of two rotating particles was achieved by 
Wan et al. [9] However, a more serious attempt must 
be directed to a suspension of finite volume fraction. It 
is therefore instructive to extend the existing theories to 
the scenario of many rotating particles, as in realistic ER 
suspensions. In which case, it is necessary to assess the 
local field effects. 

However, for an individual rotating particle in such a 
system, e.g., see Fig. 1, the relaxation of the induced 
surface charge on the particle will be towards the sum 
Etotai of the external applied field Eo and the local field 
El due to all other particles. Moreover, the presumable 
polarization P of the particles is generally not parallel to 
the applied field Eo [Fig. 1(b)]. In a heuristic example, 
since the local field is always proportional to the polar- 
ization P by virtual of the Lorentz concept, the total field 
is neither along the polarization nor along the external 
applied field in presence of any rotational motion. The re- 
sult of the competition between rotation and relaxation, 
however, must lead to the prescribed steady-state polar- 
ization. As a result, the steady-state polarization must 
be determined self-consistently and it leads to an effec- 
tive medium theory for such a system of many rotating 
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particles. Regarding the local field effects, the difficulty 
lies in the tensorial nature of the problem. More pre- 
cisely, as we mentioned, it is because the polarization is 
not parallel to the applied field due to the competition 
between the displaced surface charge on the particles via 
mechanical rotation and the polarization's relaxation to- 
wards the applied field. Therefore, a successful theory 
must accommodate these properties correctly. 

In the next section, we shall put forth an effective 
medium theory that is able to describe such properties 
of ER fluid consisting rotating suspensions in a very gen- 
eral situation. Then in Sec. Ill we examine in detail a 
special case of two-dimensional (2D) configuration and 
present some simple analytical results which will be con- 
firmed numerically. At the same time, we report an in- 
teresting off-diagonal polarization response arising from 
the dynamic effect. We conclude the paper with some 
discussions in the last section. 



II. DYNAMIC EFFECTIVE MEDIUM THEORY 

Let us consider a volume fraction p of homogeneous 
spherical particles of dielectric constant ei, embedded in 
a host medium of unit dielectric constant [see Fig. 1(a)]. 
An external electric field E is applied to the system. We 
then assume a uniform rotational angular velocity u>o for 
all the suspending particles. This could be a typical ER 
fluid under shear, [8] wherein the suspended particles are 
of micrometeric size and have a large ratio of surface-to- 
bulk polarization. 

If we first ignore the interaction between the particles, 
the dynamic equation for the induced dipole moment p 
of a single particle is given by [9] 



p + w xp 



!(p-p<°>) 



(1) 



where p(°) is the dipole moment of the particle in the ab- 
sence of rotation, and To is the relaxation time of the non- 
interacting particles. The steady-state dipole moment of 
the particle can be expressed in terms of p(°) 



,(o) 
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where uoq = \u>o\- The steady- state dipole moment p 
of Eq. (2) is generally expressed in a sum of axial and 
solenoidal terms. Furthermore, we can show that 
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In the presence of interaction, one can define the effective 
dielectric constant e in terms of the polarization P = 
X) i Pi/V of the system [20] 
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• E , 



(4) 




FIG. 1: (Color online) (a) Cartoon representation for a system 
containing rotating particles of finite volume fraction, (b) 
Schematic relationship of the single particle polarization and 
the self- consistent macroscopic polarization. 



where / is the identity tensor. Note that by definition, 
e is a tensor and frequency dependent. By invoking the 
Lorentz local electric field for the steady- state responses, 
the induced dipole moment is given by 



Pv = a 



47T 

E + y P 



(5) 



where a is polarizability tensor of the rotating particle 
which is also frequency dependent, and v is the average 
volume per particle in the suspension. By substituting 
Eq. (4) into Eq. (5), we arrive at the dynamic Maxwell- 
Garnett formula 



.I=g«.( 6 + 2Z). 



(6) 



Thus Eq. (6) represents a generalization of the static 
Maxwell-Garnett approximation (MGA) to that of ro- 
tating particles. In the absence of rotation, the polariz- 
ability tensor cx.^ = ofi'I becomes diagonal 



v(°) 



ci-1 
ei +2 



(7) 



where a is the radius of particle and e\ is the dielectric 
constant of the particles. Thus = e^I can be ob- 
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tained from the (static) Maxwell-Garnett formula: 

e(°> - 1 /ei -1\ 
7WT-2 =P {—2)> (8) 

where p = Ana 3 /3v is the volume fraction. Since p can be 
expressed in terms of p(°\ ol can be expressed in terms 
of a<°) 

fl^HMf) (0)=A (0) 
l + (^or ) 2 

I 



where u>o<^o represents a dyadic tensor. It is worthy not- 
ing that the above derivation and the resulting Eq. (9) 
are independent of coordinates and rotationally invariant 
in space. Specifically, the most general form of A reads, 
for example, in the coordinate represented by the inset 
of Fig. 1(b) 
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where uOq = 



J-'^ + o;^. Thus, the self-consistent equa- 



tion [Eq. (6)] becomes simultaneous equations for the 
components of e, which are generally quite complicated. 



III. SPECIAL CONFIGURATION IN 2D 



The solenoidal term consisting of the cross product 
uo x p(°) in the numerator of Eq. (2) indicates a dero- 
gation of the dynamic effect in the direction parallel to 
u>o- Thus Eq. (6) can be discussed in a 2D case if Eo is 
parallel to one of uo x1 uo yi and uo z . In other words, one 
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FIG. 2: (Color online) The diagonal (e xx ) and off-diagonal 
(e xy ) components of the effective dielectric constant e for the 
case of uniform rotation of angular velocity u>o = (jOqz. Note 
that there is a maximum off-diagonal response where e xx = 
exy, and the peak position ooq (marked by the vertical dashed 
lines) of e xy shifts to the lower part as the volume fraction p 
or the dielectric contrast e± increases. 



can always decompose ujq to two parts, parallel and per- 
pendicular to Eo. Due to this fact, it is instructive to 
simply consider a special geometry in which ujq = ouqz, 
while p lies in the xy plane, then Eq. (9) represents a 
transformation a = A(uoqTq) • ol^\ and 
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which indicates that the rotation does not affect the po- 
larization along z axis. In this case, the self-consistent 
equations become 
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(e xx + 2 - uoTo€ xy ), (12a) 
2 [e xy + uoT (e yy + 2)] (12b) 



c xy> ^yz 



(12c) 



0. 



c xxi £yx — ^-xyi ^ 

Thus, e zz is described by the usual MGA, as expected. 
Equations (12) suggest a complex notation. If we adopt 
complex notation p = p x + ip y , the transformation p = 
ApW can be achieved by a complex number A = 1/(1 — 
iuuoTo). If we write e = e xx + ie xy , the dynamic MGA 
[Eq. (6)] becomes a scalar but complex equation 



P A 



ei-1 
ei +2 



(13) 



One can show that Eqs. (6) and (13) lead to exactly the 
same self-consistent equations [Eqs. (12a) and (12b)] for 
this special configuration in 2D. The polarization can be 
obtained from Eq. (4) 



3/? 



= 3/3/(1-/?) 

47r(l — /3 — icuoTo) 47r(l — iuoqt) ' 



(14) 



where 



and 



= P 
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As a consequence, when we vary the rotational angu- 
lar velocity a;o, = Im[P] will exhibit a peak when 
uoqt = 1. These are clearly shown by Fig. 2, where we 
plot the diagonal effective dielectric constant e xx (solid 
lines) and the off-diagonal effective dielectric constant e xy 
(dashed lines) as functions of uuq. The effective dielectric 
constant can be solved directly from Eq. (6), Eq. (12), 
or from Eqs. (4) and (14). It is very interesting to see 
that an electric field Eo = Eqx can induce a polarization 
response in the perpendicular direction i.e., nonvan- 
ishing P y appears due to the rotational motion u>o = uoqz 
of the particles, which accumulates part of the surface 
charges, or rotates the dipole moment p, off the x direc- 
tion. Due to this transfer by rotational motion, the diag- 
onal response e xx decreases monotonically as the strength 
of the transfer effect increases (e.g., rotation speeds up), 
whereas a peak in e xy shows up. This observation of off- 
diagonal response is quite interesting, resembling a classic 
Hall effect. In Fig. 2, we also illustrate both the effects of 
the volume fraction p and the dielectric contrast e\ on the 
polarization responses, as well as on the peak shifting of 
£xy We have specifically used e\ = 2.5 in Figs. 2(a) and 
2(b) (panels at the left), while e\ = 5.0 in Figs. 2(c) and 
2(d) (panels at the right). Also we set p = 0.1 and 0.3 
for the upper panels and the lower panels, respectively. 
Combination of the effects of the dielectric contrast e\ 
and the volume fraction p is actually captured by the 
definition of f3 in Eq. (15). It is clearly seen in Fig. 2 
that increasing (3 = 0.033, 0.057, 0.1, and 0.171 leads 
to a shifting of the peak position (vertical dashed lines) 
of e xy at uol = 0.97/ro, 0.94/r , 0.90/r , and 0.83/r . 
These correspond to Figs. 2(a), 2(c), 2(b), and 2(d), re- 
spectively. From another perspective, it is noticed from 
Eq. (16) that the relaxation time will increase from non- 
interacting relaxation tim To to r due to the interactions. 
This is shown in Fig. 3, in which the solid line repre- 
sents Eq. (16) and the triangles (A) are the data of cjq 
extracted from Fig. 2. Then we can regard r as the in- 
teracting relaxation time, and the maximum off-diagonal 
response occurs when uoqt = 1, which is equivalent to 
uoqTq = 1. These results suggest a way to measure the 
interacting relaxation time in the system. 

To this end, we have demonstrated the off-diagonal 
response due to the rotational motion. We would also 
like to look upon the geometrical meaning of these re- 
sults. It is easy to tell from Eq. (14) that tan^ = c^ot, 
where <fi is the angle spanned by E and El, as shown 
in Fig. 1(b). Also by some simple calculations of the 
angles (</>, and 7) shown in Fig. 1(b) we find that 
tan# = tan(0 — 7) = tan O = c^oTo, where 6q is the 




FIG. 3: Interacting relaxation time r predicated by Eq. (16) 
compares favorably to numerical data extracted from Fig. 2. 



angle of induced dipole moment p with respect to Eo for 
the case of an isolated rotating particle at angular veloc- 
ity uoq. [9, 11] This is as expected and required by the 
self-consistency of the theory, showing that this relation- 
ship still holds for an individual rotating particle in the 
many-particle system: in Fig. 1(b) the dipole moment p 
and the total field E to tai on this particle spans the angle 
= #o- Now it is clear that the maximum off-diagonal 
response appears when <j> = 7r/4, where c^o = uJq indi- 
cating a resonance. Right at this point, the parallel and 
perpendicular effective dielectric constant are equal, i.e., 
£xx = txy, which is already demonstrated by Fig. 2. In 
the limit of uo — > 0, e xx is readily predicted by the static 
MGA [i.e., Eq. (8)] and e xy vanishes. This represents the 
case without rotation, corresponding to <p — > and — > 
in Fig. 1(b). Whereas in the opposite limit of uoq — > 00, 
both e xx and e xy approach unity (see also Fig. 2). This 
means that the rotating particle lost its tensorial polar- 
izability a of Eq. (9) [or seen from Eq. (11)] and the sys- 
tem dielectrically "sees" no polarizable suspension, i.e., 
the suspensions become "invisible" , because the very fast 
rotational motion distributes the surface charges homo- 
geneously along the circumference of the particle, accu- 
mulating no net charges at any inter facial place. In this 
case, both <j> and lost definition because of El=0. As 
a result, the system responses just as the pure host me- 
dia. It should be remarked that at a finite yet very fast 
rotation ouq^> (j) ^ tt/2 and — > 0. 

Furthermore, we examine the maximum (e.g., denoted 
by Max[- • • ]) polarization responses of the system within 
the dynamic effective medium theory, in terms of both 
diagonal e xx and off-diagonal e xy . It is quite straightfor- 
ward from the previous discussion and the data shown 
in Fig. 2 that Max^J and Max^y] occur at c^o = 
and coo = coq respectively, for a certain j3. This is also 
explicitly indicated by Eq. (14). Open circles (O) and 
squares (□) in Fig. 4 show the Max[e xai ] and Max^ J ex- 
tracted from a series figures like Fig. 2 for various /?, while 
solid lines represent e xx (cjQ = 0) and e xy (oj = (Jq). The 
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FIG. 4: The dependence of the maximum polarization re- 
sponses on f3. (a) e xx which occurs when ljo — 0, i.e., pre- 



dicted by the static MGA, and (b) e x 

CJq — Uq. 



which occurs when 



good agreements further show the usefulness of Eqs. (13) 
and (14) in the 2D case, which appears simple but con- 
tains enough information in understanding the interest- 
ing properties of the system. It is reasonable that both 
Max^J and Max^^] increase monotonically with /?, 
because commonly either an increased dielectric contrast 
or an increased volume fraction lead to larger effective 
dielectric constant. 



IV. DISCUSSION AND CONCLUSION 

To conclude, we have developed an dynamic MGA to 
take into account the many-body interactions between ro- 
tating particles in a dynamic ER fluid. The theory, like 
static MGA, stems from the application of the Lorentz 
local field concept to a system of polarized particles in 
an external electric field. The only difference is now that 
the particles are all rotating with some angular velocity 
which leads to an off-diagonal polarization. We have as- 
sumed uniform rotational angular velocity in the present 
paper. This assumption may break down in shear flow 



under non-steady conditions. If the angular velocity with 
the same magnitude uoq distributes uniformly in the solid 
angle, there is a reduction factor associated with the av- 
erage longitudinal (diagonal) response (parallel to E ) 
while the transverse (off-diagonal) responses vanish iden- 
tically. Namely, (py) = p(°) [1 - 2^r 2 /3(l +^o r o)] while 
(p_l) = 0, here (• • • ) represents average over the solid an- 
gle. 

Along with the present results, it is easy to derive a 
dynamic Bruggemann effective medium approximation 
(EM A). In contrast to the present asymmetric dynamic 
MGA, it is symmetrical and appropriate for two types 
of rotating particles (e.g., of e\ and 62) whose volume 
fractions satisfy p\+p2 = 1. In the same sprit of Eq. (13) 
in the special 2D geometry, the self-consistent equation 
of dynamic EMA reads 



Aipi 



(ci/c c + 2) 



Mp2 



feA~e ~ 1) 
(e 2 /e e + 2) 



= 0, 



(17) 



where the two complex numbers A\ and A2 correspond 
to the transformation coefficients of two types of rotat- 
ing particles. These complex numbers are assumed to 
be independent of the effective medium of e e . Physically 
Eq. (17) means that the induced dipole moment of type 
a (=1, 2) rotating particle in an effective medium points 
to a different direction due to the relaxation of surface 
charges. However, since the different types of particles 
are embedded in an effective medium, the average dipole 
moment must vanish, thus yielding the dynamic EMA 
which would be applicable to analogous problems like 
magnetotransport in granular materials. These dynamic 
effective medium theories [e.g., Eq. (6) or Eq. (17)] pro- 
vide guidelines for numerical simulation on the dynamic 
effect in ER fluids [5] and should be distinguished from 
those effective medium theories developed for compos- 
ites of anisotropic inclusions, [21] bi- anisotropic media, 
[22] and chiral media, [23, 24] etc. 

Extension to nonlinear ac responses [25, 26] is possi- 
ble. In which case, peak response can be found in parallel 
response of higher harmonics. If we further apply an ac 
field at uo along the x axis, then there are two charac- 
teristic responses at the sum and difference of frequen- 
cies: c^o + and c^o — w. When we vary uj, it can be 
shown that P y exhibits either peak at uoq + u\ or 002 — ooq. 
Thus, we can also determine the interacting relaxation 
time r by measuring the ac responses from the relations 
(c^o + oj\)t = 1 or (u)2 — uq)t = 1. It is quite interesting 
that the off-diagonal polarization shows a peak below coo 
in the linear case while it peaks exactly at uoq = 1/Vo 
for isolated rotating particle. So the ac field of uo can 
be applied perpendicular or parallel to Eo- The parallel 
component of P shows no peak (it simply decrease as uoq 
increases) in the linear case but its third harmonic term 
as a function of c^o has a peak. [25] 

Finally, we would like to discuss several complications 
that may introduce corrections to the local field and 
therefore the macroscopic responses (e.g., the effective 
dielectric constant e). (1) since most ER fluid might 
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be anisotropic, e.g., somehow chains are formed in the 
direction of applied electric field Eg, the depolarization 
term may not be isotropic. When the lattice symme- 
try is lowered by an external means — under the influence 
of an external force/torque — the lattice is deformed, ei- 
ther lengthened in one direction and/or contracted in the 
other direction. [27] These geometric anisotropy effect 
will introduce a non-isotropic local field correction. [12] 
This additional anisotropy must also be accommodate 
in a consistent way. (2) The assumption of isotropic de- 
polarization is necessary because of rotational invariance. 
More precisely, the Lorentz cavity is always spherical [see 
Fig. 1(b)]. When P is rotated away from Eo due to the 
competition between the displacement and relaxation of 
surface charges, the Lorentz relation El = 47rP/3 re- 
mains unchanged. However, the detailed dynamical cor- 
relation between the rotating particles may be further 
explored in presence of a distribution of rotational speed 
or intrinsic relaxation time. (3) There may exist certain 
uniaxial gradient in one direction in the system, [28] for 
instance, more particles aggregate in a specific part of the 
system. (4) We have not taken into account the effect of 
Brownian motion. In the presence of an external electric 



field, the Brownian motion would become weaker than 
the case of no external electric field, which is caused to 
appear by the electric polarization of the particle. How- 
ever, competition between rotational Brownian motion 
and particle fibrillation is also a concern if the particles 
are small enough to ensure rotational diffusion. 
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